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Abstract: Using two nuclear models, i) the relativistic continuum Hartree-Bogoliubov (RCHB) theory 
and ii) the Weizsácker-Skyrme (WS) model WS*, the performances of nine kinds of kernel functions in 
the kernel ridge regression (KRR) method are investigated by comparing the accuracies of describing the 
experimental nuclear charge radii and the extrapolation abilities. It is found that, except the inverse 
power kernel, other kernels can reach the same level around 0.015-0.016 fm for these two models with 
KRR method. 'The extrapolation ability for the neutron rich region of each kernel depends on the 
trainning data. Our investigation shows that the performances of the power kernel and Multiquadric 
kernel are better in the RCHB+KRR calculation, and the Gaussian kernel is better in the WS*+KRR 
calculation. In addition, the performance of different basis functions in the radial basis function method 
is also investigated for comparison. The results are similar to the KRR method. The influence of different 
kernels on the KRR reconstruct function is discussed by investigating the whole nuclear chart. At last, 
the charge radii of some specific isotopic chains have been investigated by the RCHB--KRR with power 
kernel and the WS*--KRR with Gaussian kernel. The charge radii and most of the specific features in 
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these isotopic chains can be reproduced after considering the KRR method. 
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1 Introduction 


Charge radius is one of the most basic properties 
of the atomic nucleus. With the improvement of ex- 
perimental techniques '' ^, great progress has been 
achieved in the measurement of the nuclear charge ra- 
dius, and more than 1000 experimental nuclear charge 
radii have been obtained ^H. However, this number 
is still far from enough compared to the possible exist- 
ing nuclei predicted by various nuclear theories ~’ 
Charge radius is quite important to understand vari- 
ous nuclear phenomena, such as shell structure u a 
shape phase transition Num the neutron skin and 
halo EUN etc. Therefore, accurate nuclear charge 
radius predictions are quite important. 

Usually the nuclear charge radius is described 
by the A! law with A the mass number. There- 
after, various phenomenological formulae originated 
from the A‘/? law have been proposed The 
nuclear charge radius can also be obtained by the nu- 
clear models such as the macroscopic-microscopic mod- 
els , and the microscopic models based on the 
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non-relativistic V9 and relativistic P 四 I density 
functionals. Each of these models can give a good de- 
scription of the nuclear charge radii across the whole 
nuclear chart. In addition, the local-relation-based 
models are also adopted to predict the charge 
radius of the unknown nucleus with the information of 
its existing neighbors. 

In recent years, machine learning (ML) has been 
widely used in nuclear physics . Various ML ap- 
proaches have already been adopted to improve the nu- 
clear charge radius predictions, e.g., the Bayesian neu- 
ral network , the artificial neural network isped 
the radial basis function (RBF) approach EN the k- 
ernel ridge regression (KRR) method ER EJ etc. By 
training the ML network with the charge radius resid- 
uals, i.e., deviations between the experimental and the 
calculated nuclear charge radii, the ML methods can 
reduce the corresponding root-mean-square (rms) de- 
viations to 0.01-0.02 fm. 

KRR method is one of the most popular ML meth- 
ods =I It has provided successful descriptions 
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for various aspects of nuclear physics such as nuclear 
mass , nuclear charge radius nuclear en- 
ergy density functionals , neutron-capture reaction 
cross-sections Eu etc. In these investigations, Gaus- 
sian kernel is adopted due to the good extrapolation 
power, which originates from the decay behavior of 
the Gaussian kernel with the increasing extrapolation 
distance. It can automatically identify the limit of the 
extrapolation distance and avoid the risk of overfitting 
at large extrapolation distance. 

There are a lot of commonly used kernel function- 
s in the KRR method. Recently, Wu has shown that 
the performances of the KRR approach in nuclear mass 
predictions can be affected by different kernels EJ In 
our previous work, Gaussian kernel is adopted in the 
KRR or extended KRR method in the nuclear charge 
radius predictions. So it is necessary to study the 
performance of the KRR method with different ker- 
nels in the nuclear charge radius predictions. In the 
present work, the performances of the KRR method 
for nuclear charge radius predictions with nine kind- 
s of kernel functions, including Gaussian, Laplacian, 
Matern, Cauchy, Multiquadric, inverse Multiquadric, 
Logarithm, power, and inverse power, are compared 
by training the radius residuals of relativistic contin- 
uum Hartree-Bogoliubov (RCHB) theory B and the 
Weizsücker-Skyrme (WS) model WS* V, For the ex- 
perimental charge radii of 1014 nuclei with Z > 8, the 
rms deviations for RCHB theory and WS* model are 
about 0.035 fm and 0.021 fm, respectively. The depen- 
dence of the kernel functions on the trained data will 
also be investigated. 

This paper is organized as follows. A brief intro- 
duction to the KRR method is given in Sec. 2. The 
results obtained by the KRR method with different 
kernel functions are given in Sec. 3. With different k- 
ernels, the obtained rms deviations by RCHB+KRR 
and WS*+KRR, and the extrapolation power of the 
KRR method are compared. In addition, the perfor- 
mance of different basis functions in the RBF method 
is also investigated for comparison, which is used ex- 
tensively in nuclear mass predictions with linear basis 
Ee Finally, 


function a brief summary is given in 


Sec. 4. 


2 "Theoretical framework 


The KRR function S(z;) can be written as 


m 


S(z) 2 9 K (ay, vi)w, 


i=] 


(1) 


where x; are locations of nuclei in the nuclear chart 
with zi — (Zi,Ni), m is the number of training da- 


ta, wi are weights to be determined, and K (æj, x) is 
the kernel function, which characterizes the similarity 
between the data. 
By minimizing the following loss function 
2 2 
Lw) = 3 [S(2;) - v (zT --Allell , 


i—1 


(2) 


the weights w; can be determined, where w = 
(w1,...,Wm). The hyperparameter A (A > 0) determines 
the regularization strength and is adopted to reduce 
the risk of overfitting. Minimizing Eq. (Ø) leads to 


w—(K--AI) y, (3) 


where I is the identity matrix and K is the kernel 
matrix with elements Ki; = K(xi,x;). 

In the present work, the following nine kinds of 
kernels are adopted, 


Gaussian K(r) — exp (-r?/20”) ; (4) 

Laplacian : K(r)=exp(—r/o) , (5) 

Matern : K(r)=(1+3r/5c) exp (—3r/5o(6) 

Cauchy : K(r)=1/(1+r?/o) , (7) 

Multiquadric (MQ) : K(r)= Vr? o? , (8) 

inverse MQ :  K(r)— 1/Vr2 +02 ; (9) 

Logarithm : K(r)=In(1+4+r7), (10) 

power : K(r)=r° , (11) 

inverse power K(r)=1/r° , (12) 
where the Euclidean norm 

r-|ei-zj|-v(Zi-Z;? -(Ni-N;? — (3) 


is defined to be the distance between two nuclei. The 
a (c > 0) in each kernel is a hyperparameter and it- 
s optimal value should be obtained according to the 
data of the nuclear charge radius. 

The KRR function (M) is trained to reconstruct 
the charge radius residuals AR(Z, N) = R&?(Z, N) — 
R°'(Z,N). Once the weight parameters w; are ob- 
tained, the KRR function S(Z, N) can be calculated 
for every nucleus. Therefore, the predicted charge 
radius for a nucleus (Z,N) is given by RSPR = 
R®™ (Z, N)+S(Z,N). 


3 Results and discussion 


In the present work, 1014 experimental data with 
proton number Z > 8 are considered, which are taken 
from Refs. B, E]. There are 885 data in Ref. 回 and 
129 new data in Ref. [I]. In our previous work, the 
hyperparameters A had been carefully validated using 
the leave-one-out cross-validation to be 0.02 for the 
RCHB+KRR and 0.01 for WS*--KRR P”, and these 
two values would be adopted in the present work. 


TANG Lu, ZHANG Zhenhua : Comparison of different kernel functions in nuclear charge radius predictions by the kernel 
"S XH ridge regression method 3， 


The hyperparameter o in each kernel function is 
determined by the leave-one-out cross-validation with 
given A in the KRR method. The predicted radius for 
each of the 1014 nuclei can be obtained by the KRR 
method trained on all other 1013 nuclei with given set 
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of hyperparameters c and A. The optimized hyperpa- 
rameters o (see Table 1) with given 和 can be obtained 
when the rms deviation between experimental and cal- 
culated radii gets its minimum. 
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(color online) The rms-deviations Arms between the KRR/RBF predictions with nine different 


kinds of kernel functions and the experimental data as functions of the hyperparameter c. The black 
solid lines denote the RCHB+KRR results; the red dashed lines denote the RCHB-4-RBF results; the 
blue dotted lines denote the WS*+KRR results; the olive dash dotted lines denote the WS*+RBF 


results. 


Figure 1 shows the rms-deviations Arms between the 
KRR/RBF predictions with nine different kinds of kernel 
functions and the experimental data as functions of the hy- 
perparameter c. The black solid lines denote the RCH- 
B+KRR results; the red dashed lines denote the RCH- 
B--RBF results; the blue dotted lines denote the WS*--KRR 
results; the olive dash dotted lines denote the WS*+RBEF re- 


sults. Table 1 shows the minima of rms deviations Arms and 


the corresponding optimal hyperparameters o in the leave- 
one-out cross-validation for the KRR and RBF approaches 
with nine different kernel functions as shown in Fig. 1. The 
notation “all” denotes the total 1014 data; “885” denotes 
the 885 data in Ref. B]; “129” denotes the 129 new data in 
Ref. (Ej. 

It can be seen that, except the inverse power kernel, oth- 


er kernels can reach the same level around 0.015-0.016 fm for 


these two models with both KRR and RBF methods. In fact, 
there is no improvement by the inverse power kernel in all of 
these calculations. The Gaussian kernel in the RCHB+RBF 
method is a little larger compared to other calculations. The 
results indicate that the KRR method with different kernels 
can reach similar accuracies if proper values of hyperparam- 
eters are adopted, and the results are quite similar to the 
RBF method. Therefore, in the applications of predicting 
nuclear charge radii for the nuclei that not far away from the 
experimentally known region, the prediction accuracy may 
be hardly affected by the choices of the kernel functions. We 
also can see that the predicted results for the 129 new data 
are better than the 885 data in Ref. [B]. We also can see that 
in contrary to the RBF method, the rms deviation obtained 
by the KRR method does not grow that high for very large o 
values except the MQ kernel and the Logarithm kernel. This 
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demonstrates that the penalty term in the KRR method can 
prevent overfitting even for very unrealistic kernels. 

We can see from Fig. 1 and Table 1 that in these cal- 
culations, the parameters c in MQ kernel K (r) = Vr? +0? 
are quite close to zero (especially for the WS* data set), and 


the o in power kernel K(r) = r^ are quite close to one (es- 


Tab. 1 


pecially for the RCHB data set). Both of these two kernels 
approach to the linear kernel K(r) =r with the optimized 
hyperparameter, and their results are quite close to each 
other. This means that the linear kernel adopted in KRR 
or RBF method ™ is reasonable in the investigation of the 


nuclear charge radius. 


The rms deviations Arms (in unit of fm) and the corresponding optimal hyperparameters o in the 


leave-one-out cross-validation for the KRR and RBF approaches with nine different kernel functions. 
The notation “all” denotes the total 1014 data; “885” denotes the 885 data in Ref. [B]; “129” denotes 


the 129 new data in Ref. [Hf. 


Kernels Gaussian Laplacian Matern Cauchy MQ Inverse MQ Logarithm Power Inverse Power 
RCHB+KRR 
c 2.68 4.86 2.14 1.52 0.19 0.89 1.43 0.98 - 
Arms (all) 0.0157 0.0155 0.0157 0.0154 0.0154 0.0152 0.0154 0.0155 0.0350 
Arms (885) 0.0166 0.0165 0.0167 0.0166 0.0165 0.0163 0.0164 0.0165 0.0361 
Arms (129) 0.0097 0.0102 0.0100 0.0110 0.0099 0.0105 0.0101 0.0100 0.0260 
RCHB+RBF 
c 0.86 8.16 0.51 1.48 0.20 0.89 1.45 0.99 - 
Arms(all) 0.0169 0.0155 0.0159 0.0153 . 0.0154 0.0152 0.0154 0.0155 0.0350 
Arms (885) 0.0185 0.0165 0.0171 0.0166 0.0165 0.0163 0.0164 0.0165 0.0361 
Arms (129) 0.0127 0.0101 0.0106 0.0110 0.0099 0.0105 0.0101 0.0100 0.0260 
WS*+KRR 
c 0.70 1.66 0.31 0.96 0.01 0.72 1.07 0.76 - 
Arms(all) 0.0155 0.0158 0.0157 0.0155 0.0159 0.0156 0.0158 0.0158 0.0210 
Arms(885) 0.0167 0.0169 0.0169 0.0166 0.0169 0.0166 0.0168 0.0168 0.0215 
Arms (129) 0.0141 0.0121 0.0133 0.0127 0.0114 0.0120 0.0113 0.0113 0.0177 
WS*--RBF 
c 0.70 1.62 0.31 0.95 0.01 0.72 1.08 0.76 - 
Arms(all) 0.0154 0.0158 0.0157 0.0155 0.0159 0.0156 0.0158 0.0158 0.0210 
Arms(885) 0.0167 0.0169 0.0169 0.0166 0.0169 0.0166 0.0168 0.0168 0.0215 
Arms (129) 0.0141 0.0122 0.0133 0.0127 0.0114 0.0120 0.0113 0.0113 0.0177 


To investigate the extrapolation power of the KRR and 
RBF methods with these nine different kernels, the 1014 nu- 
clei with known charge radii are redivided into one training 
set and six test sets as follows. For each isotopic chain with 
more than nine isotopes, the six most neutron-rich nuclei 
are classified into six test sets according to the distance to 
the last nucleus. Test set 6 (1) has the longest (shortest) ex- 
trapolation distance. The hyperparameters obtained by the 
leave-one-out cross-validation in Table 1 remain the same in 
the following calculations. 

Figure 2 shows the comparison of the extrapolation abil- 
ity of the KRR (left panels) and RBF (right panels) method- 
s with eight kinds of kernels and two data sets (RCHB and 


WS*) for six different extrapolation distances. The result- 


s with inverse power kernel are not shown here since there 
is no improvement by this kernel in all the calculations. It 
can be seen clearly that with these two different data sets, 
the extrapolation ability of these kernels are totally different. 
There is no overfitting in both the KRR and RBF methods 
when the radius residual in RCHB is adopted, whereas the 
overfitting appears at extrapolation distance equals to 2 for 
all the kernels in KRR and RBF methods with WS*. For 
the RCHB model, the extrapolation ability of KRR method 
with the power kernel and MQ kernel are better than oth- 
er kernels. This may be due to the charge radius residual 
in RCHB model is larger. The predicted KRR function ap- 
proaches to zero with large extrapolation distance. However, 


for the WS* model, the charge radius residual is quite smal- 
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l, overfitting appears when the extrapolation distance gets 
larger. It can be seen that, among all of these eight kernel- 
s, the results with Gaussian kernel are better both in KRR 
and RBF methods. Therefore, the present results show that 
the kernel functions or the basis functions in the KRR/RBF 
method should be checked carefully with given data set. It 
should be noted that a statistical study of the predictive pow- 


er of the KRR method is also quite important. However, un- 


RCHB+KRR 1 (b RCHB+RBF IH RCHB/WS* 
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like the Bayesian neural network, it is difficult to introduce 
the statistical methods in the KRR approach. Previous in- 
vestigations have shown that if the hyperparameters in the 
KRR method changed a little, the obtained nuclear mass 
and charge radii remain nearly unchanged. Therefore, the 


extrapolation power is not so sensitive to the slight change 


of hyperparameters. 


2 3 4 


25 PA, distance 


Fig. 2 


(color online) Comparison of the extrapolation ability of the KRR (left panels) and RBF (right 


panels) methods with eight kinds of kernels and two data sets (RCHB and WS") for six different 


extrapolation distances. 


Since the KRR method with nearly all the kernels can 
reach the same level for the charge radius prediction, it 
is quite interesting to study the influence of different k- 
ernel functions on the very neutron rich region where no 
data exists. With the optimized hyperparameters (oc,A), 
the reconstructed function S(N,Z) for every nuclei can 
be calculated by KRR method. Fig. 3 shows the differ- 
ences AR = RP 
calculated values by RCHB and WS* models (upper pan- 
els), the KRR reconstructed function S(N,Z) with power 
and Gaussian kernels (middle panels), and the differences 
AR(KRR) = R*?*? — (R° - SFRFE) between experimental and 
the predictions of these two models with KRR corrections 
3(a) and (b) that, the 
differences between experimental and the calculated values 
by RCHB model are larger than those for the WS* model. 
While after considering the KRR corrections, the two pat- 


一 R° between experimental and the 


(lower pannels). It can be seen in Fig. 


terns of AR become quite similar to each other [Fig. 3(e) 
and (f)]. Fig. 3(c) and (d) show that, the KRR reconstruct- 


ed function 9 with the power kernel remains still quite large 


value even far away from the known region. While for the 
Gaussian kernel, it become to zero for the very neutron-rich 
nuclei, which is due to the decay behavior of the Gaussian 
kernel. It means that for a given nucleus, very little infor- 
mation can be learned from the nuclei far away from it. For 
comparison, the KRR corrections considered with the Gaus- 
sian kernel for RCHB model and power kernel for WS* model 
are shown in Fig. 4(c) and (d). After considering the differ- 
ent kernels, RCHB+KRR and WS*+KRR still give quite 
similar AR patterns [Fig. 4(e) and (f)], which are also quite 
close to those in Fig. 3(e) and (f). It can be seen in Figs. 3 
and 4 that, if the Gaussian kernel is adopted, the decay be- 
havior of the KRR reconstructed function 9 is more obvious. 
'The Gaussian kernel can avoid the overfitting to a large ex- 
tent. However, sometimes the corrections provided by the 
KRR method with Gaussian kernel maybe not enough com- 
pared to other kernels (e.g., the power kernel) if the charge 
radius residuals are large. Therefore, the adopted kernel is 
essential to the predicted KRR reconstructed function S in 


the very neutron-rich region. 
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Fig. 3 (color online) The differences AR = R**P — R°* between experimental and the calculated 
values by RCHB and WS* model (upper panels), the KRR reconstructed function S(N,Z) with 
power and Gaussian kernels (middle panels), and the differences AR(KRR) = R®? — (R° + SKRR) 
between experimental and the predictions of these two models with KRR corrections (lower pannels). 
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Fig. 4 (color online) The same as Fig. 3, but the KRR corrections are considered with the Gaussian 
kernel for RCHB model and power kernel for WS* model. 
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To see the power of the KRR method in describing 
the nuclear charge radii with the favored kernel functions, 
Fig. 5 shows the comparison between the experimental and 
calculated root-mean-square nuclear charge radii for (a) Ne 
(Z — 10), (b) Ca (Z = 20), (c) Zn (Z =30), (d) Sn (Z — 50), 
(e) Gd (Z — 64) and (f) Hg (Z — 80) isotopes. The ex- 
perimental data are shown by the black solid squares. The 
calculated results by RCHB (WS*) with and without KRR 
method using power (Gaussian) kernel are denoted by blue 
(red) open and solid triangles (circles), respectively. It can 
be seen in Fig. 5 that in general, the experimental data can 
be reproduced by RCHB and WS* models. Since only spher- 
ical shape is considered in the RCHB model, the deviation 
between the data and the calculated results is larger than 
the WS* model for these six isotopes. After considering the 
KRR method, the results of these two models are both im- 
proved and they are quite similar to each other. It can be 


seen that the odd-even staggering missing in these two mod- 


3.05 3.60 


m Exp. 
人 -RCHB -O- WS* 
3.55 H-A- RCHB + power KRR 4.00 
—e- WS* + Gaussian KRR 


els appears with the KRR corrections, especially for the Ca 
[Fig. 5(b)] and Hg [Fig. 5(f)] isotopes. It should be noted 
that in most cases, the phases of the calculated odd-even 
staggerings are opposite to the data and the amplitude is 
also not quite well reproduced. Our previous investigation 
shows that after including an odd-even term in the kernel 
function in the KRR method, the experimental OES could 
be reproduced quite well EJ Tt also can be seen in Fig. 5(d) 
that in the RCHB model, the calculated charge radii change 
smoothly when across the neutron N = 82 shell, and it can 
not reproduce the abrupt kink at N — 82. After the KRR 
method being considered, this kink can be reproduced quite 
well by the RCHB model. For the proton rich side of Ne iso- 
topes, the charge radii increase drastically with decreasing 
'This is 
Ex ux 


neutron number, especially for the "Ne [Fig. 5(a)] . 
due to the existence of the proton halo in this nuclei 
It can not be reproduced by the KRR method, which is due 


to the insufficient of the training data for the halo nuclei. 
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(color online) Comparison between the experimental and calculated root-mean-square nuclear 


charge radii for (a) Ne (Z = 10), (b) Ca (Z —20), (c) Zn (Z —30), (d) Sn (Z =50), (e) Gd (Z = 64) 


and (f) Hg (Z — 80) isotopes. 


The experimental data are shown by the black solid squares. The 


calculated results by RCHB (WS*) with and without KRR method with power (Gaussian) kernel 
are denoted by blue (red) open and solid triangles (circles), respectively. 


4 Summary 


In this work, the performances of nine kinds of kernel 
functions in the KRR method are investigated by compar- 
ing the accuracies of describing the experimental nuclear 
charge radii and the extrapolation abilities. Two data sets 
are adopted, i.e., the RCHB theory and the WS* model. It 


is found that, except the inverse power kernel, other kernels 


can reach the same level around 0.015-0.016 fm for these two 
models with KRR method both with the charge residuals in 
RCHB theory or WS* model. The results are quite similar to 
the RBF method. The results indicate that the KRR/RBF 
method with different kernels can reach similar accuracies 
It has 


been found that the extrapolation ability for the neutron 


if proper values of hyperparameters are adopted. 
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rich region of each kernel depends on the trained data. Our 
investigation shows that the performances of the power ker- 
nel and Multiquadric kernel are better in the RCHB+KRR 
and RCHB--RBF calculations. 
rameters, these two kernels are approximate to the linear 
kernel K(r) =r. In the WS*+KRR and WS*+RBF calcu- 
At last, the charge 
radii of Ne, Ca, Zn, Sn, Gd and Hg isotopic chains have 
been investigated by the RCHB+KRR with power kernel 
and the WS*+KRR with Gaussian kernel. The charge radii 
and most of the specific features in these isotopic chains can 


With the optimal hyperpa- 


lations, the Gaussian kernel is better. 


be reproduced after considering the KRR method. 
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- 10. 原子 核 物 理 评 论 OX 


核 岭 回归 方法 中 采用 不 同 核 函 数 对 于 电荷 半径 预测 的 比较 


EAM, 张振华 1,21 
( ! 华 北 电力 大 学 数理 学 院 ， 北 京 102206 


?河北 省 物理 学 与 能 源 技 术 重点 实验 室 ， 华 北 电力 大 学 ， 保 定 , 071000 ) 


摘要 : ”通过 采用 i) 相对 论 连 续 谱 Hartree-Bogoliubov (RCHB) 理论 以 及 ii) Weizsücker-Skyrme (WS) #2! WS", 
o d mc ac ss B. 
本 工作 发 现 ， 在 使 用 KRR 方法 对 于 电荷 半径 的 描述 时 ， 除 了 逆 圭 律 核 ， 采 用 其 他 核 函 数 计 算得 到 的 电荷 半径 的 均 
La ucc uc li ai 
表明 在 RCHB-KRR 的 计算 中 ， 震 律 核 以 及 多 元 二 次 核 的 外 推 效 果 较 好 ; 在 而 WS'-KRR 的 计算 中 ， 高 斯 核 的 
外 推 效 果 较 好 。 此 外 ， 本 文 也 研究 了 在 径 向 基 函 数 方法 中 采用 不 同类 型 的 基 函 数 对 于 电荷 半径 的 描述 能 力 ， 其 结果 
与 KRR 方法 类 似 。 随 后 ， 通 过 对 整个 核 素 图 进行 计算 ， 研 究 了 不 同 核 函 数 对 于 KRR 重 构 函 数 的 影响 。 最 后 ， 利 用 
考虑 竹 律 核 的 RCHB 十 KRR 以 及 高 斯 核 的 WS* 十 KRR 方 法 对 一 些 特定 同位 素 链 中 原子 核 的 电荷 半径 进行 了 研究 。 考 
虑 KRR 方 法 后 ， 能 够 很 好 地 再 现 这 些 同位 素 链 中 原子 核 的 电荷 半径 以 及 其 中 的 一 些 特征 。 
关键 词 : 电荷 半径 ; 机 器 学 习 ; 核 岭 回归 


Wis AHA: XXXX-XX-XX 
基金 项 目 : 国家 自然 科学 基金 资助 项 目 (11875027) 
+ 通信 作者 : IKIRE, E-mail: zhzhang@ncepu.edu.cn 


